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Abstract. The properties of starbursts in ultraluminous infrared galaxies are dis- 
cussed, with particular emphasis on the fueling, the amount of extinction and the 
intrinsic properties of the nuclear starbursts. It is shown by the example of NGC 6240 
that the H2 vibrational lines can be used to measure the rate of gas inflow into the 
potential well, which is sufficient to fuel a nuclear starburst of the intensity required 
to account for the far-infrared emission. It is shown that in Arp220 the faintness of all 
tracers of ionized gas can be accounted for by Lyman continuum absorption by dust 
within the ionized regions, combined with significant (but not extreme) extinction; 
there is no reason to invoke the presence of extreme extinction, an old starburst, or an 
additional non-stellar power source in Arp220. 



1 Introduction 

Are ultraluminous infrared galaxies (ULIGs) powered by intense bursts of star 
formation or is an additional source of energy, such as an active galactic nucleus 
(AGN) required? In order to address this question, it is instructive to compare 
the near-infrared spectra of ULIGs and lower luminosity starburst galaxies. In 
starburst galaxies of low or moderate luminosity, the dominant emission line is 
the Br7 line, underlining the importance of massive young stars in the energetics 
of these objects. The brightest H2 rovibrational line, the H 2 v = 1—tO S(l) line, 
is typically fainter @,|,|^JllJ||Q . In contrast, in ULIGs the H 2 v = l->0 



S(l) line is significantly brighter than the Br7 line, which is often not even 
detected WM- This behaviour shows that ULIGs are not just scaled- up starburst 
galaxies. An extreme case of is presented by NGC 6240 (Fig. @), where the H 2 
v — 1— >0 S(l) line is 40 times brighter than Br7. The low ratio of Br7 luminosity 
to FIR luminosity in ULIGs has been used to argue against star formation as 
the power source of ULIGs, for instance in the nearby ULIG Arp220, where a 
starburst with the low Br7 luminosity observed can account for at most 10% of 
the bolometric luminosity of the galaxy pp2| . However, spectroscopy of ULIGs 
at longer wavelengths with the Infrared Space Observatory (ISO) revealed bright 
emission lines from powerful starbursts, which are the dominant power source 
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Fig. 1. Spectra of NGC6240 in selected parts of the near- infrared H and K bands, 
integrated over a 4'.'4 region j^] 

in most of the objects studied, but are obscured at shorter wavelengths Q. 
Motivated by these results, in this paper the physical processes revealed by the 
i\-band spectra of ULIGs are reexamined in the context of the starburst scenario. 
In Sect. H the origin of the H2 line emission is discussed, while Sect. || addresses 
the faintncss of Br7. 

2 Origin of the H 2 emission in ULIGs 

The most extreme vibrational H2 emission is found in the nearby luminous 
merger NGC6240: 7 • 10 7 L Q is emitted in the H 2 v = l->0 S(l) line alone 
(for Hq = 75 km s" 1 Mpc -1 and with no correction for extinction). This line 
contains 0.012% of the bolometric luminosity of NGC 6240, which is considerably 
higher than any other galaxy [fLs) . Together the vibrational lines may account 
for 0.1% of the total bolometric luminosity. 

Imaging of the H2 v — 1— >0 S(l) emission from NGC 6240 has shown that 
the H2 emission peaks between the two remnant nuclei of the merging system 
[ fl8| . This morphology provides a unique constraint on the excititation mecha- 
nism, since it argues against any scenario where the excitation is dominated by 
the stellar component (e.g., UV-pumping, excitation by shocks or X-rays from 
supernova remnants). Instead, the favoured excitation mechanism is slow shocks 
in the nuclear gas component, which, as shown by high resolution interferometry 
in the CO J = 2->l line (]]|, also peaks between the nuclei of NGC 6240. 
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Fig. 2. High resolution imaging of NGC6240 in the H 2 v = l-»0 S(l) line and the 
2.15 Lim continuum with NICMOS/HST p[ 



What is the role of these shocks? In the shocks mechanical energy is dissipated 
and radiated away, mostly in spectral lines (principally H2, CO, H2O and [Oi] 
lines). This energy is radiated away at the expense of the orbital energy of the 
molecular clouds in the central potential well. Consequently, the dissipation of 
mechanical energy by the shocks will give rise to an infall of molecular gas to the 
centre of the potential well. Therefore, the H2 vibrational lines measure the rate 
of infall of molecular gas into the central potential well. This conclusion can be 
quantified by writing 



ad 



(i) 



where L ra< a is the total luminosity radiated by the shocks and Ldis the dissipation 
rate of mechanical energy, giving rise to a molecular gas infall rate Mh 2 given 

by 



1 • , 

Ldis = ^Mb_ 2 v 



(2) 



where v is the circular orbital velocity at the position where the shock occurs. 

Using a -fsT-band extinction of 0™15 [p| , the total luminosity of H 2 vibrational 
lines from NGC 6240 becomes 7.2 • 10 8 L Q ; inclusion of the purely rotational 
lines observed with ISO approximately doubles this number, so that £ ra d = 
1.5 • 10 9 L Q . 

In order to use this number to estimate Mh 2 , it is necessary to establish more 
accurately the fraction of the H2 emission that is due to infalling gas. Observa- 
tions with NICMOS on the Hubble Space Telescope (HST) provide the required 
information (Fig. ^) p9[ . The NICMOS image shows that the emission consists 
of a number of tails (presumably related to the superwind also observed in Ha 
emission), and concentrations assocated with the two nuclei, and a further con- 
centration approximately (but not precisely) between the two nuclei. The relative 
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brightness of the H2 emission from the southern nucleus is deceptive, since this 
nucleus is much better centred in the filter that was used for these observations 
than the other emission components, in particular the northern nucleus. Taking 
this effect into account, it is found that 32% of the total H2 flux is associated 
with the southern nucleus, 16% is associated with the northern nucleus, and 
12% with the component between the two nuclei, the remaining 40% being asso- 
ciated with extended emission. Using inclination-corrected circular velocities of 
270 and 360 km s _1 for the southern and northern nucleus respectively Jl^j , and 
of 280 km s~ x for the central component the mass infall rates derived using 
are 80 M Q yr _1 for the southern nucleus, 22 M© yr _1 for the northern 
nucleus and 28 M Q yr _1 for the central component. 

The derived molecular gas inflow rate to the two nuclei is remarkably close 
to the mass consumption rate by star formation of approximately 60 M yr , 
indicating that the H2 emission from the nuclei directly measures the fueling of 
the starbursts in these regions. This analysis shows that the central regions of 
NGC 6240 are being fueled at a rate sufficient to maintain starburst activity at 
the level required to account for the FIR luminosity. 

3 Dusty, compact starbursts in ULIGs 

Accepting the starburst model, the faintness of Br7 and other recombination 
lines remains to be addressed. As pointed out in Sect, [l], it is evident that selective 
extinction towards the regions of most recent massive star formation plays a 
significant role in suppressing the Br7 emission. Can this effect be quantified? 

The nearby ULIG Arp 220 has been studied in detail with the ISO satelite. 
Based on an upper limit on the ratio of the [S ill] 18 and 33 iim lines, and a 
high Bra over Br7 flux ratio, an extinction Ay = 50 ± 10 m , located purely in an 
absorbing foreground screen was proposed H]. However, the supporting argu- 
ments have now weakened significantly. In the first place, a better understanding 
of the calibration of the [S ill] 33 |4.m spectrum has made the upper limit on the 
[S III] 18 to 33 (J,m flux ratio less strict by approximately a factor of two Sec- 
ondly, the Bra flux from ISO |Q is almost certainly an overestimate: the Bra 
line displays, on top of a higly structured baseline, a double-peaked structure, 
which is absent in any other line (including long-wavelength lines such as the 
well-detected [S in] 33 11m line). The velocity difference between the two peaks in 
Bra is approximately 600 km s _1 and therefore cannot be attributed to motion 
of the two nuclei of Arp 220, which have a radial velocity difference of approxi- 
mately 200 km s _1 The same velocity difference of 200 km s" 1 is found in 
long-slit Br7 spectra of the Arp 220 nuclei j|]. An integrated high-resolution Br7 
spectrum (Fig. ||) shows no trace of a double-peaked structure, indicating that 
the high Bra flux found with ISO is most likely dominated by structure in the 
spectral baseline. The extinction derived from the Bra/Br7 ratio should thus 
be used as an upper limit. A further argument against an obscuring foreground 
screen with Ay ~ 50 m is furnished by the derived Lyman continuum fluxes, 
which increase as as shorter wavelength tracers are used 0, a behaviour sug- 
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Fig. 3. Integrated spectrum of Br7 emission from Arp220 ]l7| 

gesting that the extinction has been overestimated. Finally, a strong limit on the 
presence of an intense, highly obscured, but otherwise normal starburst follows 
from the upper limit to the free-free emission at millimetre wavelengths 
These results indicate that either a foreground screen with a lower visual extinc- 
tion, or a model with mixed emission and absorption needs to be adopted. 

This result does however not imply the presence of a strongly aged starburst, 
or an additional source of power in Arp 220, since the preceding analysis ignored 
the role of Lyman continuum absorption by dust within the ionized regions. 
If most of the ionizing radiation is absorbed by dust grains rather than hy- 
drogen atoms, a dust-bounded (rather than hydrogen-bounded) nebula results, 
and all tracers of ionized gas (recombination lines, fine-structure lines, free-free 
emission) will be suppressed. If the H n regions in Arp 220 are principally dust- 
bounded, the observational properties of Arp 220 can be accounted for, even with 
only moderate extinction. Since the dust would also absorb far-ultraviolet radi- 
ation longwards of the Lyman limit, the formation of photon-dominated regions 
would also be suppressed, and the thus the same mechanism can account for the 
faintness of the 158 |im [C n] line in Arp 220 and other ULIGs ||. 

Is the Arp 220 starburst dominated by dust-bounded H n regions? The av- 
erage molecular gas density in the ~ 10 10 M nuclear molecular complex in 
Arp 220 is n^ 2 ~ 2 • 10 4 cm~ 3 The strong emission from high dipole mo- 
ment molecules such as CS, HCO + and HCN argues for even higher densities: 
~ 10 10 Mq of molecular gas (i.e., all of the gas in the nuclear complex) has a 
density 7Jh 2 ~ 10 5 cm -3 fD| , At such densities the ionized nebulae created by 
hot stars are compact or ultracompact Hii regions, where 50 to 99% of the Ly- 
man continuum is absorbed by dust |2^ | . Observationally, hydrogen-bounded and 
dust-bounded Hn regions can be distinguished by the quantity R = L-pm/ Lr^'- 
for a wide range of parameters, R < 3570 implies that the nebula is hydrogen- 
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bounded, while R > 35700 implies that the nebula is dust-bounded B. For 
Arp220, the Br7 luminosity of 1.3 • 10 6 L Q (from the spectrum in Fig. H with 
a distance of 77 Mpc) implies R = 1.6 • 10 5 assuming an obscuring foreground 
screen with Ay — 20 m (a model consistent with the results discussed above). 
Even with a foreground extinction of Ay — 50 m (which is most likely an over- 
estimate, as discussed above), a ratio R = 1 • 10 4 would result, so that even in 
that case the absorption of Lyman continuum radiation by dust would play a 
significant role. The star formation takes place in (ultra)compact Hil regions, 
where all of the usual tracers of ionized gas (recombination lines, fine-structure 
lines, free-free emission) are quenched, not extincted. While this result signifi- 
cantly complicates the interpretation of diagnostics of massive star formation in 
ULIGs, it is save to conclude that the properties of Arp 220 can be accounted for 
by an intense, and significantly (but not extremely) obscured starburst. There is 
no reason to invoke the presence of extreme extinction, a strongly aged starburst, 
or an additional power source in Arp 220. 
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